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ABSTRACT 



We present cosmological simulations which closely mimic the real Universe within ^ 100 Mpc 
of the Local Group. The simulations, called Constrained Simulations, reproduce the large-scale 
density field with major nearby structures, including the Local Group, the Coma and Virgo 
2 ! clusters, the Great Attractor, the Perseus-Pices, and the Local Supercluster, in approximately 

*— > ' correct locations. The MARK III survey of peculiar velocities of the observed structures inside 

jrt . 80/i~^Mpc sphere is used to constrain the initial conditions. Fourier modes on scales larger 

then a (Gaussian) resolution of w i?g = 5ft.~^Mpc are dominated by the constraints, while 
small scale waves are essentially random. The main aim of this paper is the structure of the 
^^ ' Local Supercluster region (LSC; '^ 30/i~^Mpc around the Virgo cluster) and the Local Group 

H \ environment. We find that at the current epoch most of the mass (« 7.5x 10^^/i~^Mq) of the LSC 

is located in a filament roughly centered on the Virgo cluster and extending over ~ 40/i~"'^Mpc. 
The simulated Local Group (LG) is located in an adjacent smaller filament, which is not a part 
of the main body of the LSC, and has a peculiar velocity of ~ 250 km s~^ toward the Virgo 
cluster. The peculiar velocity field in the LSC region is complicated and is drastically different 
from the field assumed in the Virgocentric infall models. We find that the peculiar velocity flow 
in the vicinity of the LG in the simulation is relatively "cold" : the peculiar line-of-sight velocity 
dispersion within 7/i~^Mpc of the LG is < 60 km s~^, comparable to the observed velocity 
dispersion of nearby galaxies. 

Subject headings: cosniology:theory - large-scale structure - methods: numerical 

1. Introduction testing theoretical models in the highly nonlinear 

regime. However, comparisons of simulations with 

Over the last two decades cosmological simu- observational data are typically done only in a sta- 

lations have proved to be an invaluable tool in ^^^^^^^1 g^^g^ The standard approach is to assume 

; a cosmological model and to use the appropriate 

Hubble Fellow , „ , , . , . , ; i , • 

power spectrum of the primordiaf perturbations 



to construct a random realization of the density 
field within a given simulation volume. The evo- 
lution of the initial density field is then followed 
using a numerical code and results are compared 
with observations. A wide variety of statistical 
measures has been used for comparisons, including 
the two-point correlation function, the power spec- 
trum, the mass function, different shape statistics, 
etc. 

The statistical approach works well if there is a 
statistically representative sample of objects with 
well understood selection effects for both the ob- 
served Universe and the simulations. The number 
of sample objects should be sufficient to overcome 
the cosmic variance. The approach fails, however, 
in cases dealing with rare objects such as the Great 
Attractor or specific configurations as, for exam- 
ple, in the problem of the local peculiar velocity 
field (Peebles 1992). The traditional solution has 
been to choose "observers" or regions in the simu- 
lations that resemble the desired configuration as 
closely as possible. The results of numerical sim- 
ulations aimed at identifying structures similar to 
those observed locally in random realizations of 
the initial density field were generally inconclu- 
sive because one is never sure that the selection of 
objects in simulations is right. 

The catalog of galaxy groups constructed using 
the CfA redshift survey (Davis et al. 1982) is a 
good illustration of the difficulties one encounters 
in comparing model predictions to observations 
(Huchra & Geller 1982; Nolthenius et al. 1997). 
While the data-set was sufficiently large for statis- 
tical comparisons with models, a large fraction of 
groups in the sample is in one large object: in the 
Local Superclustcr (LSC). For detailed compar- 
ison with observational data it was important to 
identify groups in a similar environment in cosmo- 
logical simulations. For example, it was crucial to 
mimic the environment of the Local Group (LG) 
in the LSC (e.g., Nolthenius et al. 1997). This was 
done by simply choosing an "observer" 20ft-^^Mpc 
away from a cluster with mass comparable to that 
of the Virgo Cluster. Unfortunately, it is not clear 
what "similar environment" actually means and 
that simply placing the "observer" at some dis- 
tance from "Virgo" cluster resolves the issue. 

Another problem that is not easy to address 
using statistical comparisons is the peculiar ve- 
locity field around the LG. The observed low ve- 



locities of nearby galaxies have long being per- 
ceived as one of the important challenges for struc- 
ture formation models (Peebles 1992). Because we 
deal with rather specific environment, in simula- 
tions we must find objects, which are analogous 
to the real Local Group. The interpretation of the 
galaxy velocity field is further complicated by the 
possibility of velocity bias (Gelb & Bertschinger 
1994; Colin et al. 2000). Nevertheless, it was sug- 
gested that, even after taking the uncertainties 
into account, observed peculiar velocities of galax- 
ies around the LG are unusually low if compared 
to the typical velocities of galaxies in cosmological 
simulations (Schlegel et al. 1994; Governato et al. 
1997). 

There are two possible ways of addressing prob- 
lems, which require comparison with specific as- 
tronomical objects or with specific environments. 
(i) We may simply ignore them and wait for larger 
observational data-sets, such as the Sloan Digital 
Sky Survey, and pursue the statistical comparison 
of the observed data with the theoretical predic- 
tions. Clearly, the increasing amount and quality 
of data will continue to fuel the statistical way 
of testing the cosmological models. At the same 
time, observational data always will be more ex- 
tensive and more complete for nearby, somewhat 
unique objects. Thus, by doing only statistical 
tests we are bound to lose information on the best 
observed astronomical objects in our sample. 

(ii) An alternative approach, adopted in the 
present study, is to find a way of making sim- 
ulations, which reproduce the large scale struc- 
tures of the real Universe, while keeping the ini- 
tial conditions consistent with the power spectrum 
of a given cosmological model. In particular, we 
want to simulate the small scale nearby structures 
within the correct large environment of the real 
Universe as revealed by various large scale sur- 
veys. This is done by setting initial conditions 
of the simulations by constrained realizations of 
Gaussian fields (Hoffman & Ribak 1991), thereby 
constructing the density and velocity fields which 
agree both with the observed large scale structures 
and with the assumed theoretical model. This 
approach is called here as constrained simulation 
(CS). 

Low resolution CSs were made before by Kolatt 
et al. (1996) and Bistolas & Hoffman (1998), us- 
ing the IRAS 1.2Jy redshift survey to set the con- 



straints, with the main aim of producing a semi- 
reahstic non-Hnear reaHzation of the large scale 
structure. Recently, such CSs were made using the 
MARK III survey of radial velocities (Willick et al. 
1997), and were used (van de Weygaert & Hoff- 
man 1999, 2000) to study the origin of the local 
cold velocity field. The study presented here ex- 
tends this work by using numerical simulations of 
much higher spatial and mass dynamic range. The 
dynamic range achieved by the Adaptive Mesh Re- 
finement simulations used in this study far exceeds 
the range of all previous CSs. The combination 
of CSs and adaptive mesh refinement is the opti- 
mal strategy for simulating formation of the local 
structures at high spatial and mass resolution in 
a computational box sufficiently large to properly 
model the large-scale tidal field. 

The present work is the first in a series of pa- 
pers studying the dynamics of the nearby universe 
by means of CSs. The main purpose of this pa- 
per is to present the general approach and details 
of the method. The specific CS analyzed here is 
a relatively low resolution one, while forthcoming 
CSs will use higher resolution initial conditions 
and will reach higher mass and spatial dynamic 
range. As an illustration, this paper addresses the 
question of morphology of the density field and 
dynamics in the LSC region. In the following sec- 
tions we present analysis of the structure of the 
Local Supercluster, as well as the smoothness of 
the Hubble expansion and peculiar velocity field in 
the simulated LSC region. We will also discuss the 
accuracy of the Virgoccntric model in describing 
the velocity field in our simulation. The upcoming 
studies will focus on a variety of problems, includ- 
ing gas-dynamic simulations of the Virgo cluster 
and the LSC region (Kravtsov et al. 2001), high- 
resolution simulation of the LG, halo/galaxy for- 
mation in the Local Void region, and biased galaxy 
formation in the local universe. 

The paper is organized as follows. A brief 
review of studies of the LSC region is given in 
§ 2. The formalism of constrained realizations, 
the MARK III survey, and the construction of the 
initial conditions are presented in § 3.1, § 3.2, and 
§ 3.3, respectively. The details of the numerical 
simulation are described in § 3.4 and § 4. The 
analysis of the CS is presented in §5 and we con- 
clude the paper with a discussion of the results in 
§6. 



2. The Local Supercluster 

Ever since the first large samples of nebulae 
compiled by William and John Herschel, it was 
known that there is a marked excess of bright ob- 
jects in the northern hemisphere around the Virgo 
cluster, with most nebulae concentrated in a wide 
band spanning some 100° (e.g., Reynolds 1923; 
Lundmark 1927; Shapley 1934)^. In 1950s de Vau- 
couleurs (1953, 1958) was the first to argue that 
this excess corresponds to a real 3-dimensional 
structure in galaxy distribution. At first de Vau- 
couleurs (1953) called the structure the "Super- 
galaxy" (and introduced the supergalactic coordi- 
nates in analogy with galactic coordinate system) 
conjecturing that the apparent flattened galaxy 
distribution represented a large-scale galaxy-like 
system, later he subsequently changed the name 
to the "Local Supercluster" (LSC; de Vaucouleurs 
1958), which has been used in the literature ever 
since. 

The debate on whether the LSC is a physi- 
cal system or a chance alignment of galaxies has 
continued over several decades (e.g., Bahcall & 
Joss 1976; de Vaucouleurs 1976, 1978) and was 
settled with the advent of large redshift surveys 
of nearby galaxies (Sandage & Tammann 1981; 
Fisher & TuUy 1981; Davis et al. 1982). Yahil 
et al. (1980) and TuUy (1982), for example, an- 
alyzed the morphology of the 3-dimensional dis- 
tribution of galaxies in the Sandage & Tammann 
(1981) and Fisher & TuUy (1981) surveys, respec- 
tively, and showed convincingly that the concen- 
tration of galaxies along the supergalactic plane 
on the sky does indeed correspond to a flattened 
large-scale structure in the distribution of nearby 
galaxies. Einasto et al. (1984) applied a battery of 
tools to quantify the statistical properties of the 
galaxy distribution in the LSC region and com- 
pared the results with quantitative and qualita- 
tive predictions of several cosmological models. In 
particular, they showed that the LSC is similar 
in structure and morphology to other known su- 
perclusters such as the Perscus-Pisccs supercluster 
and is a part of the large-scale network of clusters, 
sheets, and filaments. 

More recent studies (TuUy & Fisher 1987; 



^See interesting historical note on "discovery" of the LSC 
by de Vaucouleurs (1989) and § 2 in historical review by 
Biviano (2001) 



Karachentsev & Makarov 1996; Lahav et al. 
2000) showed that the main body of the LSC 
is a filamentary structure that extends over some 
40/i~^Mpc and is roughly centered on the Virgo 
cluster. The whole region is dominated by sev- 
eral clusters (Virgo, Ursa Major, and Fornax are 
the most prominent), groups, and filaments, the 
latter bordering nearby voids (such as the Local 
Void). The LG is located in the outskirts of this 
region in a small filament extending from the For- 
nax cluster to the Virgo cluster (Karachentsev & 
Makarov 1996). The LSC contains a fair number 
of galaxy groups (Huchra & Geller 1982) and its 
skeletal structure is well traced by radio galaxies 
and AGNs (Shaver & Pierre 1989; Shaver 1991). 

The quest for the Hubble constant and the ad- 
vent of peculiar velocity surveys of galaxies fueled 
studies of galaxy dynamics in the LSC region (see, 
e.g., Davis & Peebles 1983; Huchra 1988, and ref- 
erences therein). In particular, it was recognized 
that the infall of the LG toward the LSC could 
be used to measure the mass-to-light ratio (and 
hence the matter density in the Universe, Oq) on 
large scales (Silk 1974; Peebles 1976). The Vir- 
gocentric infall model (Gunn 1978) was used to 
estimate ^Iq ^ 0.2 — 0.3 from the observed local pe- 
culiar flow of galaxies (Huchra 1988). This model 
requires rather large non-linear corrections (Yahil 
1985; Villumsen & Davis 1986) and was shown to 
be quite inaccurate when tested using outputs of 
numerical cosmological simulations (Villumsen & 
Davis 1986; Lee et al. 1986; Cen 1994; Governato 
et al. 1997), rendering the conclusions about Og 
inconclusive. 

The quiescence of the local peculiar velocity 
field (de Vaucouleurs 1958; Sandage & Tammann 
1975; Rivolo & Yahil 1981; Sandage 1986; Gi- 
raud 1986; Schlegel et al. 1994; Karachentsev & 
Makarov 1996; Sandage 1999; Ekholm et al. 2001) 
is a long standing puzzle which presents a chal- 
lenge for models of structure formation (e.g., Pee- 
bles 1992). That this puzzle is by no means solved 
is clear from the recent study by Sandage (1999) 
who comments that "the explanation of why the 
local expansion field is so noiseless remains a mys- 
tery." The radial peculiar velocity dispersion of 
galaxies within w 5/i~^Mpc around the Milky Way 
is only ^ 50 — 60 km s~^ and the local Hubble 
flow agrees with the global expansion law on large 
scales to better than 10% (Sandage 1999). A strik- 



ing manifestation of the coldness of the flow is the 
fact that outside the Local Group there are no 
nearby galaxies with blueshifts with respect to the 
Milky Way. 

Clearly, successful cosmological models should 
provide a plausible explanation for such a "cold" 
velocity field. However, we are dealing with a 
unique and limited set of data which makes it 
difficult to employ the usual statistical compar- 
isons between models and observations. Theo- 
retical studies employing cosmological simulations 
concluded that the CDM models fail to natu- 
rally reproduce the "coldness" of the local veloc- 
ity field (Schlegel et al. 1994; Governato et al. 
1997). Schlegel et al. (1994) used a variety of 
criteria to identify LG candidates in their sim- 
ulations of the standard CDM and mixed DM 
(MDM, or cold+hot dark matter) models. They 
found that the observed "coldness" of the local 
fiow cannot be reproduced in the then standard 
CDM model, but could be plausibly reproduced 
in the MDM model. Governato et al. (1997) used 
a set of somewhat different criteria to identify 
LG counterparts in the simulations of the stan- 
dard and open CDM models. They selected bi- 
nary systems consisting of halos with circular ve- 
locities of 125 — 270 km s^^ approaching each 
other, which do not have other massive neighbors 
within 3 Mpc. An additional subsample was con- 
structed by requiring that the selected binary sys- 
tems are located at 5 — 12/i^^Mpc from a Virgo 
sized cluster. This study confirmed that the stan- 
dard CDM model cannot reproduce the observed 
low value of velocity dispersion. The velocity dis- 
persions around groups in the open CDM model 
with mass density Oq = 0.3 (~ 150 — 300 km s^^) 
were found to be lower than those in the SCDM 
(~ 300-700 km s^i), but still higher than the ob- 
served value. The authors concluded that "neither 
the r^ = 1 (CDM) nor n = 0.3 (OCDM) cold dark 
models can produce a single candidate LG that 
is embedded in a region with such small peculiar 
velocities." 

The conclusions of these studies, however, have 
to be viewed with caution as they explored only 
a limited range of cosmological models (in partic- 
ular, the currently favored ACDM model was not 
considered) and have used a statistical approach to 
select possible LG candidates from a fairly small- 
size simulation volume. For example, it is not clear 



whether simple criteria lead to selection of objects 
representative of the LG or, if the LG environment 
is rare, that the box size was sufficient to contain 
a suitable counterpart. It is also not clear what is 
the largest scale that influences the specific prop- 
erties of the LG environment. The local "cold" 
flow may be influenced by the mass distribution 
on much larger scales and be induced by coupling 
of the small and large scales via the velocity shear 
(see van de Weygaert & Hoffman 1999, 2000). 
The local criteria for the definition of LG counter- 
parts are then clearly inadequate. 

3. Constrained Simulations 

3.1. Wiener Filter and Constrained Real- 
izations of Gaussian Fields 

The standard cosmogonical framework of the 
CDM family of models assumes that structures 
evolve out of small perturbations in an expand- 
ing Friedmann universe. It is assumed that these 
perturbations constitute a Gaussian random field. 
Redshift and radial velocity surveys provide in- 
formation that enables the reconstruction of the 
large-scale structures (LSS) in the density field of 
the nearby Universe. An efficient algorithm for 
reconstructing the density and velocity fields from 
sparse and noisy observations of the LSS, such as 
the one provided by redshift and velocity surveys, 
is provided by the formalism of the Wiener filter 
(WF; for a general overview see Zaroubi et al. 
1995). The application of the WF requires some 
model for the power spectrum that defines the 
statistical properties of the perturbation field. A 
Bayesian approach to the problem of choosing the 
"best" model consists of finding the most probable 
model over an assumed parameter and/or model 
space by means of maximum likelihood analysis, 
and using this model for calculating the WF (see, 
however, Hoffman & Zaroubi 2000 for a discussion 
of the limitations of the Bayesian approach.) The 
application of the general formalism to the case of 
radial velocity surveys follows Zaroubi et al. (1997) 
and Zaroubi et al. (1999). 

Below we briefly describe the WF reconstruc- 
tion method. Consider a survey of N objects with 
radial velocities {Mi}i=i....,Ar, where 

Ui = v(rj) -ri +ei, (1) 

Here v is the three dimensional velocity, r^ is the 



position of the i-th data point, and ei is the statis- 
tical error associated with the i-th radial velocity. 
Given that the assumed model, and that the sta- 
tistical errors are well understood, the data auto- 
covariance matrix can be readily evaluated: 



UzUjJ = fj{|v(r,)v(rj)^fj + , 



(2) 



Here 



denotes an ensemble average. The last 

term erf, is the error covariance matrix. The ve- 
locity covariance tensor is calculated using linear 
theory. 

Given a peculiar velocity dataset, the WF pro- 
vides the minimum variance estimation of the un- 
derlying field(s). In the case where the underly- 
ing field is Gaussian, the WF-reconstructed den- 
sity field coincides with the most probable and the 
mean fields for the data, and also with the maxi- 
mum entropy estimation (Zaroubi et al. 1995). As- 
suming the linear theory relation of the radial ve- 
locity with the full three-dimensional velocity and 
density fields, the WF estimation of these fields is: 



WF 



(r) = (v{r)u,j(u,ujj 



and 



6^'ir)^(6{r)u,)(u,u 



(3) 



(4) 



The (v{r)ui) and (S{r)ui) are the cross radial 

velocity ~ three-dimensional velocity and density 
correlation matrix. 

The WF is a very conservative estimator. In the 
absence of "good" data, namely where the data 
is sparse and/or noisy, it attenuates the estimate 
towards its unbiased mean field, which in the cos- 
mological case is the null field. Thus, by construc- 
tion the WF suppresses some of the power that is 
otherwise predicted by the assumed model. The 
WF often produces an estimated field that is much 
smoother than the typical random realization of 
the assumed power spectrum would be. In partic- 
ular, the WF estimator is not statistically homo- 
geneous. A way of providing the missing power 
and regaining the statistical homogeneity consis- 
tent with the data and with the theoretical model 
is provided by the method of constrained realiza- 
tions of Gaussian fields (Bertschinger 1987; Hoff- 
man & Ribak 1991, 1992). The constrained re- 
alizations provide a realization of the underlying 



field made of two components. One is dictated 
by the data and by the model and the other is 
random in such a way that in places, where the 
WF suppresses the signal, the random component 
compensates for it. 

The WF attenuation can be overcome by the 
unbiased minimum variance estimator of Zaroubi 
(2001), where some of the missing power is pro- 
vided by the observational errors. It provides an 
attractive tool for recovering the LSS, but it is not 
suitable for providing realizations that are consis- 
tent with a theoretical model. 

The Hoffman & Ribak (1991) algorithm of con- 
strained realizations provides a very efficient way 
of creating typical realizations of the residual from 
the WF mean field. The method is based on creat- 
ing random realizations of the density and veloc- 
ity fields, (5(r) and v(r), given an assumed power 
spectrum and a proper set of random errors e^. 
The random realization is then "observed" just 
like the actual data to yield a mock velocity data 
set Ui. Constrained realizations of the dynamical 
fields are then obtained by 



OR/ 



(r) =v(r) + (v(r)wA(u,Mj) (uj-Mj) (5) 



and 

6''\r)^d{r) + (^S{r)u,)(^u,u,y\u,-u,). (6) 

The variance of the constrained realizations 
around the WF mean field provides a measure of 
the amount by which they are constrained by the 
data. This variance can be calculated rigorously 
using the auto- and cross-correlation matrices de- 
fined here (Zaroubi et al. 1995). However this cal- 
culation becomes impractical for large grids as it 
involves the inversion of large matrices. A much 
simpler and more efficient way is to construct an 
ensemble of constrained realizations and to calcu- 
late its scatter (Zaroubi et al. 1999). 

The auto- and cross-covariance matrices in the 
above equations are computed within the frame- 
work of the linear theory. The WF and the con- 
strained realizations are performed assuming that 
both the sampled data and the evaluated fields are 
in the linear regime. Indeed, for velocity data pro- 
cessed by the grouping algorithm, the linear the- 
ory provides a good approximation (see Kudlicki 



et al. 2001). The choice of the resolution at which 
the WF and fields of constrained realizations are 
evaluated is arbitrary and it can be controlled by 
a Gaussian smoothing of radius Rg and of course 
by the grid scale over which the fields are con- 
structed. Besides the practical limitations im- 
posed on choosing the grid size, there is one ba- 
sic consideration that dictates the choice of the 
grid and the smoothing. This can be best un- 
derstood in terms of the Fourier space presenta- 
tion, where in the case of unconstrained realiza- 
tions the Fourier modes are statistically indepen- 
dent. The imposed constraints introduce mode 
coupling, where the range of the coupling depends 
on the kind of data (say, density or velocity), its 
quality (errors, sparseness and space coverage) and 
the nature of the power spectrum (Bertschinger 
1987; Hoffman & Ribak 1991, 1992). In a cos- 
mological framework, only waves in a finite range 
of wavenumbers are constrained, while longer and 
shorter waves are unaffected by the data. In the 
WF reconstruction this would mean that these 
waves are set to zero, while for the constrained re- 
alizations these waves are effectively sampled out 
of unconstrained realizations. Thus, ideally the 
resolution is set to R„ 



"^T^^rnax^ where fc„ 



cor- 



responds to the maximum k constrained by the 
data. 

The WF and constrained realizations are con- 
structed assuming that the linear theory is valid 
on all scales. Thus, in principle it can be done 
with any desired resolution, in particular on scales 
that at present lie deep in the non-linear regime. 
In other words, the WF and constrained realiza- 
tions provide a reconstruction or a realization of 
how the present-day structure would appear if the 
linear theory had been valid. This limitation can 
be turned into an advantage and used as a tool 
for recovering the initial conditions that seeded 
the growth of structures in the nearby universe. 
The basic idea of the present paper is to use the 
data on scales where linear theory is applicable 
to recover the large scale fluctuations and to sup- 
plement them with fluctuations due to a random 
realization of a specific power spectrum on small 
scales. These fluctuations are extrapolated back 
in time using the linear theory to provide a recon- 
struction of the initial conditions. 



3.2. Observational Data and Prior Model 

The MARK III catalog (Willick et al. 1997) 
has been compiled from several data sets of spiral 
and elliptical/so galaxies with the direct TuUy- 
Fisher and the £)„ — a distances. The sample con- 
sists of ~ 3400 galaxies and provides radial veloc- 
ities and inferred distances with fractional errors 
~ 17 — 21%. The sampling covers the whole sky 
outside ±(20 — 30)° of the Zone of Avoidance. It 
has an anisotropic and non-uniform density that is 
a strong function of distance. The sampling of the 
density field is generally good to about 60/i~^Mpc, 
although this limit changes from 40/i~'^Mpc to 
80/i~^Mpc for some directions. The data are cor- 
rected for the Malmquist biases. A grouping pro- 
cedure has been applied to the data in order to 
lower the inhomogeneous Malmquist bias before 
it is used and to avoid strong non-linear effects, 
in particular in clusters of galaxies. This yields a 
dataset of distances, radial peculiar velocities, and 
errors for k, 1200 objects, ranging from individual 
field galaxies to rich clusters. 

The cosmological model assumed here is the 
currently popular flat low-density cosmological 
model (ACDM) with f^o = 1 -^A = 0.3, where VIq 
is the cosmological density parameter of the non- 
relativistic matter and Oa measures the cosmolog- 
ical constant A in units of the critical density. The 
Hubble constant is assumed to be /i = 0.7 (mea- 
sured in units of 100 kms^^Mpc^^) and the power 
spectrum is normalized by erg — 0.9. The ACDM 
model assumed here is consistent with all current 
observational constraints (e.g., Wang et al. 2001). 
This model is also consistent with the radial veloc- 
ity surveys including the MARK III, although the 
data favor a slightly higher value of fio (Zaroubi 
et al. 1997). 

A detailed analysis of the LSS reconstructed 
from the MARK III survey was presented by 
Zaroubi et al. (1999). The most robust features 
of the structure recovered from the MARK III 
are the Great Attractor (GA), the Perseus-Pisces 
(PP) supercluster, the filamentary LSC connect- 
ing GA and PP, and the Local Void. The nearby 
structures can be resolved using current methods 
down to a resolution of about 5/i^^Mpc. 



3.3. Reconstruction using the MARK III 
Survey 

The WF/Constrained Realizations algorithm 
has been applied to the MARK III database as- 
suming a flat ACDM model (see previous sec- 
tion). This was done on a 128"^ grid with an 
(1.25/i^^Mpc)'^ cell, thus reconstructing the den- 
sity and the velocity field within a box 160/i~^Mpc 
on a side centered on the LG. The WF is evaluated 
by calculating the various correlation functions as 
a one-dimensional k integral over the power spec- 
trum, and therefore the WF is obtained without 
imposing any boundary conditions on the grid. 
On the other hand, the constrained realizations 
are based on an unconstrained realization gener- 
ated by an FFT algorithm, which imposes peri- 
odic boundary conditions. The density and ve- 
locity fields reconstructed using the WF method 
are presented in Fig. 1, where the overdensity, 
5 = 5p/ p^ and the velocity, v, fields are evaluated 
with the Gaussian smoothing of Rg — 5ft.^^Mpc. 
A constrained realization with the same resolution 
is shown in Figure 2. 

Visual comparison of Figures 1 and 2 reveals 
that there are some similarities and some differ- 
ences. On large scales the velocity fields are more 
indicative: they show that the same large struc- 
tures are present in both distributions. Thus, the 
Great Attractor and a large void below it, and the 
Perseus-Pisces Supercluster are clearly seen in the 
velocity field (the PP supercluster does not appear 
in the density field because in this realization the 
supercluster is outside the shown slice). The Lo- 
cal Supercluster is present in both plots as a small 
"island" slightly above and to the left of the cen- 
ter. There is a small extension towards negative 
SGY, which in the real Universe hosts the Fornax 
poor cluster and the Eridanus cloud. 

Nevertheless, there are some visible differences 
between Figures 1 and 2. The WF field produces a 
smooth regular field, by attenuating a large frac- 
tion of the power predicted by the power spec- 
trum. This attenuation increases with the dis- 
tance from the LG, corresponding to noisier and 
sparser data. As a result, small-scale structures 
are found only where the data are sufficiently ac- 
curate, as, for instance, in the GA region (around 
[SGX.SGY] « [-40,0]/i-iMpc. This explains 
why in the peripheral parts of Figure 2 there are 
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Fig. 1. — Wiener filter reconstruction of tfie density 
and velocity fields from the MARK III survey of ra- 
dial velocities, evaluated with a 5h~ Mpc Gaussian 
smoothing. The plots show the structure of the su- 
pergalactic plane. The density field {top panel) is pre- 
sented by contours of constant overdensity with spac- 
ings of (5 = 0.25 and the velocity field {bottom panel) 
is shown by streamlines. The LG is in the middle 
at [0,0]. The two largest structures are the Perseus- 
Pisces Supercluster on the low right corner and the 
Great Attractor is in the middle left of the plots. 
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Fig. 2. — A constrained realization of the density and 
velocity fields constructed using constraints from the 
MARK III survey and assuming ACDM model. The 
presentation follows Fig. 1, except here the velocity 
field is not smoothed. In this realization the Perseus- 
Pisces supercluster is located slightly below the super- 
galactic plane (at SGZ « — 10/i^^Mpc) and thus is 
missing in the density plot, but it is apparent in the 
velocity field. Due to periodic boundary conditions the 
realization is missing a bulk motion of ~ 300 km s^^ 
in the direction of the Shapley Concentration, which is 
present in the WF velocity field in Figure 1. The distri- 
bution has an overall displacement of about 12/i~^Mpc 
to positive SGY with respect to the WF field. 



many more small structures. The two velocity 
fields differ in one significant aspect: the con- 
strained realization field clearly reflects the peri- 
odic boundary conditions while WF field does not. 
The large-scale bulk velocity is therefore missing 
in the case. WF velocity field has a relatively large 
bulk flow. This is why the flow around the Perseus 
cluster in constrained realization is more symmet- 
ric (infall velocities from all directions) and the 
flow around the real Perseus cluster is asymmet- 
ric: it has net velocity in the direction of the Great 
Attractor. Another effect of the periodic condi- 
tions is the shift of the whole matter distribution 
in the direction of positive SGY by 11.5/i~^Mpc 
and by — 4/i~^Mpc in SGX direction, compared 
to the WF fleld. This is especially obvious if 
one compares the position of the Virgo cluster 
in WF and in the constrained realization data. 
In order to compensate for this, we displace the 
whole distribution in supergalactic coordinates by 
Ar= [3.5, -11.5, 0.9]/i-iMpc. 

One conclusion that follows from this analysis 
is that the constrained realization adds the small 
scale power, while retaining the large scale struc- 
tures dictated by the data and manifested by the 
WF. Generally, the position and amplitude of ob- 
jects on scales of a few Mpc do vary with the dif- 
ferent constrained realizations but the large scale 
(> 10/i~"'^Mpc) structures are reproduced robustly. 

To study the constraining power of the data 
within the framework of the assumed ACDM 
power spectrum an ensemble of ten constrained 
realizations has been generated. The constrained 
realizations have been constructed with Gaussian 
smoothing of _Rg = 5 and 10ft-~^Mpc, and with 
no smoothing at all. We then estimate the mean 
variance {{Scr — Swf)'^) of the constrained re- 
alizations (the scatter around the WF 6 field) 
for all regions inside a sphere of radius R from 
the LG. By comparing the smoothed density field 
of a constrained realization to the original WF 
field, we can gauge the power of the observational 
constraints in the realization. Figure 3 shows 
the variance normalized by the unconstrained ex- 
pected variance \<5acdm)- The plot shows that 
down to a resolution of a few Mpc the region 
within « 20/i^^Mpc, which is the subject of the 
present study, is well constrained by the data: the 
variance is small with only 20-25% of the power 
contributed by the random component of the con- 
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Fig. 3. — The mean variance {{Scb. — Swf)^) of the 
ensemble of constrained realizations of the density field 
with respect to the density field obtained using WF 
reconstruction as a function of depth for a Gaussian 
smoothing of Rg = 10/i~^Mpc (dotted fine), 5/i~^Mpc 
(solid line), and no smoothing at all (dashed line). The 
variance is normalized by the variance of the ACDM 
unconstrained field smoothed with the same filter. 



strained realizations. This ratio increases with 
the depth; at large distances the realizations be- 
come unconstrained. On the other hand, on the 
scales of < 5 — 10/i^^Mpc, the fiuctuations are 
heavily dominated by the random realization and 
are virtually unaffected by the constraints (the 
field is unconstrained on scales below the grid cell, 
< 1/i^^Mpc). It follows that the unsmoothed con- 
strained realization provides us with a realization 
that has all the power demanded by the assumed 
model, with the large scales being constrained by 
the data, while the small-scale waves are largely 
unconstrained. 

3.4. Constrained Simulations 

The aim of this paper is to perform A'^-body 
simulations that match the observed local uni- 
verse as well as possible. Namely, we are inter- 
ested in reproducing the observed structures: the 
Virgo cluster, the Local Supercluster and the LG, 
in the approximately correct locations and embed- 
ded within the observed large-scale configuration 



dominated by the GA and PP. constrained real- 
izations, however, show some sizeable deviations 
from the WF-reconstructed field and the degree 
with which a given realization matches the ob- 
served density and velocity fields varies. These 
deviations are, of course, much smaller than in the 
case of an unconstrained realization, but we would 
nevertheless like to concentrate the computational 
effort on a realization that matches the observed 
configuration most closely. 

To this end, we have generated five constrained 
realizations and selected a realization (shown in 
Fig. 2) that provides the best match to the ob- 
served density field reconstruction. Specifically, 
we have performed low-resolution simulations of 
each constrained realization. These simulations 
were run with 128"^ equal mass particles with a 
particle mass of 1.6 x 10^^/i~^Mq and the high- 
est formal spatial resolution of 20/i^^kpc. All 
the simulations have reproduced the Virgo cluster, 
the Great Attractor and the Perseus-Pices super- 
cluster, with the Virgo cluster being embedded in 
a LSC. The position of the Virgo cluster varied 
in the simulations along the LSC by more than 
10/i~^Mpc either in the direction of the GA or to 
the PP supercluster. We have selected the real- 
ization in which the Virgo cluster at the present 
epoch was in an approximately correct location. 
This was the only selection criterion as the reso- 
lution of the simulations was barely sufficient to 
resolve the Virgo cluster. The small-scale struc- 
ture within the LSC region was not resolved. 

Using the selected constrained realization, two 
simulations with increasing force and mass resolu- 
tion in the region around the Virgo Cluster were 
performed. The initial conditions for these sim- 
ulations were set using multiple mass resolution. 
At the z = output of the low-resolution run, we 
selected all particles within a sphere of 25/i^^Mpc 
radius centered on the Virgo cluster. The mass 
resolution in the lagrangian region occupied by 
the selected particles was increased and additional 
small-scale waves from the initial ACDM power 
spectrum of perturbations were added appropri- 
ately (see Klypin et al. 2001, for details of the 
method). For the two high-resolution simulations, 
the particle mass in the LSC region is 8 and 64 
times smaller than in the low-resolution simula- 
tion. The highest resolution simulation has parti- 
cle mass 2.51 x IQ^h^^Mp, and the maximum for- 



mal force resolution was 2.44/i^^kpc in the LSC 
region. The results of both high-resolution simu- 
lations agree well with each other at all resolved 
scales. Below we present results only from the 
highest resolution run. 

4. Numerical simulations 

The ART TV-body code (Kravtsov et al. 1997; 
Kravtsov 1999) was used to run the numerical sim- 
ulation analyzed in this paper. The code starts 
with a uniform grid, which covers the whole com- 
putational box. This grid defines the lowest (ze- 
roth) level of resolution of the simulation. The 
standard Particles-Mesh algorithms are used to 
compute the density and gravitational potential 
on the zcroth-level mesh. The code then reaches 
high force resolution by refining all high density 
regions using an automated refinement algorithm. 
The refinements are recursive: the refined regions 
can also be refined, each subsequent refinement 
having half of the previous level's cell size. This 
creates a hierarchy of refinement meshes of differ- 
ent resolution, size, and geometry covering regions 
of interest. Because each individual cubic cell can 
be refined, the shape of the refinement mesh can 
be arbitrary and effectively match the geometry 
of the region of interest. This algorithm is well 
suited for simulations of a selected region within 
a large computational box, as in the constrained 
simulations presented below. 

The criterion for refinement is the local density 
of particles: if the number of particles in a mesh 
cell (as estimated by the Cloud-In-Cell method) 
exceeds the level nthresh, the cell is split ("refined") 
into 8 cells of the next refinement level. The re- 
finement threshold may depend on the refinement 
level. The code uses the expansion parameter a 
as the time variable. During the integration, spa- 
tial refinement is accompanied by temporal refine- 
ment. Namely, each level of refinement, I, is in- 
tegrated with its own time step Aa; = Aao/2', 
where Aao is the global time step of the zeroth 
refinement level. This variable time stepping is 
very important for accuracy of the results. As the 
force resolution increases, more steps are needed 
to integrate the trajectories accurately. Extensive 
tests of the code and comparisons with other nu- 
merical A'^-body codes can be found in Kravtsov 
(1999) and Knebe et al. (2000). 
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The current version of the ART code has the 
abihty to handle particles of different masses. In 
the present analysis this ability was used to in- 
crease the mass (and correspondingly the force) 
resolution inside a region centered around the 
Virgo cluster. The multiple mass resolution is im- 
plemented in the following way. We first set up a 
realization of the initial spectrum of perturbations 
in such a way that initial conditions for a large 
number (1024^) of particles can be generated in 
the simulation box. Coordinates and velocities of 
all the particles are then calculated using all waves 
ranging from the fundamental mode k = lix jL to 
the Nyquist frequency fc = 2-k jL x N'^^^ /2, where 
L is the box size and N is the number of parti- 
cles in the simulation. Some of the particles are 
then merged into particles of larger mass and this 
process can be repeated for merged particles. The 
larger mass (merged) particle is assigned a veloc- 
ity and displacement equal to the average velocity 
and displacement of the smaller-mass particles. 

The simulations presented here were run using 
256"^ zeroth-level grid in a computational box of 
160/i~^Mpc. The threshold for cell refinement (see 
above) was low on the zeroth level: rithrosh(O) = 2. 
Thus, every zeroth-level cell containing two or 
more particles was refined. This was done to pre- 
serve all small-scale perturbations present in the 
initial spectrum of perturbations. The thresh- 
old was higher on deeper levels of refinement: 
^thresh — 3 and rithrosh = 4 for the first level and 
higher levels, respectively. 

For the low resolution runs the step in the 
expansion parameter was chosen to be Aao = 
3 X 10""^ on the zeroth level of resolution. It was 
Aao = 2 X 10~^ for the high resolution runs. This 
gives about 500 steps for particles located in the 
zeroth level for an entire run to z = and 128,000 
for particles at the highest level of resolution. 

5. Results 

5.1. Morphology of matter distribution 

Figure 4 shows the density and velocity fields 
in a slice centered on SGZ = 0. The fields are 
smoothed with the Gaussian filter of 5/i~^Mpc 
smoothing length. No additional projection was 
used. The plane of the slice was rotated by 
30 degrees around the vertical line SGX = 0. 
This slice goes through both the Great Attrac- 



tor and the Perseus-Pisces supercluster. All ma- 
jor structures (the LSC, Great Attractor, Perseus- 
Pisces supercluster, and Coma cluster) observed 
within 100/i~^Mpc around the Milky Way exist 
in the simulations. The positions and morphol- 
ogy of these structures is, of course, fairly well 
constrained by the constraints imposed on the ini- 
tial conditions. In this projection the LSC is an 
elongated structure some 15/i~^Mpc above the ori- 
gin. It extends over ~ 40/i~^Mpc along the SGX 
axis. There is a low-density "bridge" (of overden- 
sity just above the average density), which con- 
nects the LSC with the Perseus-Pisces Superclus- 
ter. There is a also an even weaker filament con- 
necting the LSC with the Great Attractor (it is 
not obvious in Fig. 4 but is visible on the left side 
of Fig. 5). 

Figure 5 shows a zoom- in view of 45/i^^Mpc re- 
gion around the Virgo cluster. It is clear that the 
matter distribution is very far from being spher- 
ically symmetric. The main body of the Local 
Supercluster is a filament roughly centered on the 
Virgo Cluster and extending from the Ursa Ma- 
jor cluster to a concentration of several massive 
groups (region around (-10,10)/i^^Mpc). Smaller 
filaments connect the LSC to other nearby struc- 
tures. Two filaments in the upper half of the plot 
(SGX« -16/i^iMpc and « lAh-'^U-pc) connect 
the LSC with the Great Wall and the Coma clus- 
ter. The filament extending diagonally from the 
Ursa Major down to the bottom right corner of 
the slice forms a "bridge" to the Perseus-Pices 
Supercluster. The weak filament extending al- 
most straight down along the SGY axis from the 
Virgo cluster connects it to the LG and, eventu- 
ally, to the Great Attractor. The large structure 
at the bottom of the figure at {SGX, SGY) = 
(— 8, — 10)ft.~^Mpc is the location of the Fornax 
cluster counterpart in the simulation. It was just 
outside of the region of high mass and force res- 
olution. The sharp decline in the number density 
of particles in the bottom left and right corners 
of the Figure 5 marks the boundary of the high 
resolution region. 

Just as in the real Universe, the LG is located in 
a weak filament extending between the Virgo and 
Fornax clusters. This filament is a counterpart of 
the Coma-Sculptor "cloud" in the distribution of 
nearby galaxies (see TuUy & Fisher 1987, for the 
survey of nearby structures) . Figure 5 shows that 
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Fig. 4. — Density and velocity fields in a slice through the central regions of the simulated volume {SGZ — 0). The 
fields are smoothed with the 5/i~^Mpc Gaussian filter. Positions of several known objects are marked (GA - Great 
Attractor, LS - Local Supercluster, PP - Perseus- Pisces Supercluster, Coma - Coma cluster). The thick large arrow 
in the top right corner corresponds to the velocity of 1000 km s~^. The thin density contours correspond to 1, 2, 3, 
and 4 times of the average density of the matter in the Universe, while the thick contours shows the overdensity 5. 
Our position is at the origin of the system of coordinates at (0,0). Due to a rather large smoothing length, most of 
the structures prominent in this figure are well constrained to match the observed structures in the Mark III catalog. 
An important difference between the simulated and observed velocity fields is that the tidal field produced by the 
matter outside the simulated volume was removed. 
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our Galaxy is not located in a common galaxy en- 
vironment. A typical galaxy is very likely to be 
found in the main body of the LSC few mega- 
parsecs away from a large group or a cluster. At 
the same time, our location is not very special. 
There are other small filaments in LSC, which are 
similar to our Coma-Sculptor "cloud" . 

Figure 6 shows a zoom-in view of the immedi- 
ate environment of the simulated LG. Note that 
the structures at these scales are only weakly af- 
fected by constraints imposed on the initial con- 
ditions. Several possible counterparts to existing 
objects (e.g., the MW and M31, M51, NGC253) 
are marked, but their existence is largely fortu- 
itous. As can be seen in this figure, the simulated 
LG is located in a rather weak filament extending 
to the Virgo cluster (see Fig. 5). This filament 
borders an underdense region visible in the right 
lower corner of Figure 6, which corresponds to the 
Local Void in the observed distribution of nearby 
galaxies. Note that the velocity field around the 
LG is rather quiet. The peculiar velocity field in 
the Local Void exhibits a uniform expansion of 
matter out of this underdense region, while veloc- 
ities between the LG and the Virgo (upper half of 
Fig. 6) show a coherent flow onto the main body 
of the LSC. 

As can be seen from Figures 4-6 most of the 
mass in the LSC region is concentrated in a fil- 
ament centered on the Virgo cluster, while the 
LG is located in an inconspicuous albeit slightly 
over dense region neighboring the Local Void. The 
matter distribution is thus not spherically sym- 
metric as is often assumed in the Virgocentric ve- 
locity fiow models and the mass within the Virgo 
cluster itself constitutes only a small fraction of 
the LSC's mass. In the next section we will dis- 
cuss the properties of the peculiar velocity field in 
this region and compare it to the observed velocity 
of nearby galaxies. 

5.2. Peculiar velocity field in the LSC re- 
gion 

Figures 4-6 show that the peculiar velocity field 
in the LSC region is quite complex. The deviations 
from the uniform Hubble flow can be thought of as 
two different components: the bulk flows (coher- 
ent large-scale fiows which can be seen in Figs. 4-6) 
and the small-scale velocities in and around col- 
lapsed (or collapsing) objects. In underdense re- 



gions the bulk velocities exhibit a roughly spheri- 
cally symmetrical pattern typical of an expanding 
void (Bertschinger 1985): the peculiar velocities 
are steadily increasing from the center of the void 
to its bordering filaments. In Figure 5 this pattern 
can be seen clearly in the two voids in the top (be- 
tween the Virgo and the Coma clusters) and the 
lower right corner (the Local Void). Note that 
in the vicinity of the high-density regions peculiar 
velocities tend to have direction perpendicular to 
the nearest pancake or filament (see, for example, 
velocity field around the filament near the "Ursa 
Major cluster" in the middle right of Fig. 5). This 
behavior can be well understood in terms of the 
Zel'dovich (1970) pancake solution which predicts 
velocity flow in the direction perpendicular to the 
largest dimension of the collapsing pancake. 

It is also clear that peculiar velocities do not 
exhibit a simple Virgocentric spherically symmet- 
ric infall pattern usually assumed in the models of 
the local velocity field (e.g., Aaronson et al. 1982; 
Huchra 1988). Only in the immediate vicinity of 
the Virgo cluster (within ^ 3 Mpc of the clus- 
ter center) the velocity infall is Virgocentric and 
is close to being spherically symmetric. On larger 
scales the velocity fiow is roughly perpendicular to 
the filament which constitutes the main body of 
the LSC, reflecting the fact that most of the mass 
is distributed in this filament rather than being 
concentrated in one cluster. 

The predictions of the Virgocentric infall model 
fail to describe the properties of the simulated pe- 
culiar velocity field accurately both qualitatively 
and quantitatively. Figure 7 shows spherically av- 
eraged density, radial velocity, and velocity disper- 
sion profiles centered on the Virgo cluster"^. The 
average density contrast at the distance of the LG 
in the simulation is 5p/ p ~ 0.69. The dot-dashed 
line in Figure 7 shows the density model used by 
Huchra (1988) in the Virgocentric infall model: a 
constant density with added r"^ law, which is nor- 
malized to have an overdensity of 4 at the LG dis- 
tance. This model overpredicts the mass within 
the MW distance around the Virgo cluster in our 
simulation by a factor of three. This is consistent 



•^The simulated Virgo cluster has the virial radius 
960/i^^kpc indicated by short vertical dashed line in Fig- 
ure 7, the virial mass of M^-„ = 1.03 X lO^'^/i^lMo, and 
the concentration of C = 8.7 (see Navarro ct al. 1997, for 
definition). 
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Fig. 5. — Density and velocity fields in the 45/i~^Mpc region around the Virgo cluster. The velocities are in the 
Virgo cluster rest frame. The fields were smoothed with the Gaussian filter of lAh^^Mpc smoothing length. The 
circle at the origin of the coordinates marks the position of the Milky Way galaxy. Points show dark matter particles 
(10% of all particles is shown) in a slice of 10/i^^Mpc thickness centered on the Supergalactic plane {SGZ = 0). 
Contours show the projected density in the slice: the thick contour corresponds to the average matter density of the 
Universe. The thin contours mark overdensities 2, 4, 6, and so on. The length of the thick arrow in the top right 
corner shows a velocity of 500 km s^^. 
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Fig. 6. — Density (contours) and velocity (arrows) fields (smoothed with a Gaussian filter of 0.7/i^^Mpc smoothing 
length) around the LG. The slice shown has a size and thickness of 15/i^^Mpc and 5/i~^Mpc, respectively, and is 
centered on the supergalactic plane {SGZ — 0) . Points show positions of the dark matter particles (all particles are 
shown). Contours show the projected density in the slice: the thick contour corresponds to the overdensity 3; the 
other contours mark overdensities 1 and 2. The length of the thick arrow in the top right corner corresponds to a 
velocity of 500 km s^^. The velocities are plotted in the Virgo cluster rest frame. Although the matter distribution 
bears some resemblance to the observed distribution of galaxies, observational data constraints on these scales are 
weak. Note the empty region to the right of the LG. 
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with the results of Villumsen & Davis (1986), Cen 
(1994), and Governato et ah (1997) who concluded 
that the Virgocentric infall model is quite inaccu- 
rate when tested using numerical simulations. 

Figure 7 also shows that the average radial ve- 
locity of the dark matter is near zero within the 
virial radius of the Virgo cluster, which indicates 
that cluster matter is in virial equilibrium. The 
deviations of the average velocity from zero at 
larger radii shows that the matter near and just 
outside the virial radius is not fully relaxed. At 
even larger radii the profile shows the infall of un- 
virialized matter and can be well approximated 
as 145(13/i"^Mpc/r)i/2 km s"^ (shown by the 
dashed curve). The radial velocity dispersion pro- 
file in the top panel of Figure 7 is approximately 
flat at large radii (indicating the velocity disper- 
sion typical in the LSC region), while within the 
virial radius the velocity dispersion profile has a 
shape typical for the halos formed in the CDM 
cosmologies (Navarro et al. 1997; Lokas & Mamon 
2001; Klypin et al. 1999). 

Figure 8 shows the Hubble diagram for the sim- 
ulated DM halos in the Milky Way restframe. Sim- 
ilarly to observations (e.g., Sandage 1999), galactic 
halos follow the global Hubble expansion remark- 
ably well. Note that deviations from the Hubble 
flow in the vicinity of the LG (r < 7/i^^Mpc) are 
< 100 km s^^ and the flow is rather "cold". In 
particular, similarly to the observed local Hub- 
ble flow (e.g., Giraud 1986; Ekholm et al. 2001) 
only one halo (besides the M31 counterpart) has 
negative radial velocity (i.e., is blueshifted). Pe- 
culiar velocities increase at larger (r > 7/i^^Mpc) 
radii due to the large concentration of mass in the 
main body of the LSC, in the immediate vicinity 
of the Virgo cluster. Table 1 presents the radial 
and 3-dimensional velocity dispersions in differ- 
ent radial shells around the MW for halo samples 
constructed using various selection criteria. The 
columns show the shell radius (1), radial velocity 
dispersion (2), 3-dimensional velocity dispersion 
(3), number of halos in the radial shell, halo sam- 
ple selection criteria (5, with Np and Vdrc denoting 
the number of particles per halo and the circular 
velocity of the halo, respectively). The rms radial 
velocity at r = 5ft.~^Mpc is about 60 km s^^ (vir- 
tually independent of the sample selection criteria) 
and increases to « 110 km s^^ at r = 10/i^^Mpc. 

To illustrate deviations from spherical symme- 
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Fig. 7. — Spherically averaged profiles of density {bot- 
tom panel), average radial velocity {top panel, lower 
curve), and radial velocity dispersion {top panel, up- 
per curve) constructed around the center of the Virgo 
cluster. Bottom panel: The density profile is shown in 
units of the average density of the universe, po. The 
short vertical dashed line on the x-axis at --^ l/i~^Mpc 
indicates the virial radius of the Virgo Cluster. The 
dashed curve shows the NEW profile (Navarro ct al. 
1997) with an added constant po to account for the 
fiattening of the profile at large radii. The dot-dashed 
curve shows the density model assumed by Huchra 
(1988) to model the velocity field in the LSC re- 
gion. Top panel: note that the average radial veloc- 
ity of the dark matter is very small within the virial 
radius of the Virgo cluster, which indicates an es- 
tablished virial equilibrium within the cluster. The 
dashed curve outside the virial radius shows a fit to 
the infall velocity profile outside the Virgo cluster: 
145(13ft"^Mpc/r)^/^ km s"\ 



try in the velocity flow around the Virgo clus- 
ter. Figure 9 shows the peculiar radial veloci- 
ties of DM particles along the line-of-sight pass- 
ing through the the Virgo Cluster. All DM par- 
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Fig. 8. — The Hubble diagram (radial velocity as 
a function of radius) for hales in the simulation in 
the LG restframe. Individual halos are denoted by 
open circles with the circle radius representing the 
maximum circular velocity of the halo (larger circles 
correspond to more massive halos). The thick line 
shows the Hubble law, hh = -ffo?", for the value H^ = 
70 km s^^Mpc^^ assumed in the simulation, and the 
thin lines correspond to velocities «h ± 100 km s~^. 
Note that deviations from the Hubble flow in the vicin- 
ity of the LG (r < 7/i^^Mpc) are < 100 km s"^ (i.e., 
the flow is "cold"). 



tides within 1.5ft.~^Mpc around the line-of-sight 
are shown and particle velocities are computed 
in the restframe of the Virgo cluster. The most 
prominent feature in the radial velocity profile is 
the "finger of god" effect within the LG and the 
Virgo cluster. This well-known effect is due to 
large random velocities of particles within these 
massive halos. The LG and the Virgo cluster are 
the only massive systems along this line-of-sight. 
Therefore, the radial velocity dispersion of par- 
ticles outside the virialized regions measure the 
"temperature" of the intergalactic velocity field. 
Note that the typical deviations from the coher- 
ent infall are very small on both sides of the Virgo 
Gluster: - 50 - 60 km s~^ at r < 7/i~^Mpc and 
< 10 km s^^ at r > 18/i~"'^Mpc, which reflects 
the relative "coldness" of the flow. The flow is es- 
pecially cold at large radii because the matter at 



these radii undergoes a fairly uniform expansion 
out of a void. 

The solid curve in Figure 9 represents a model 
for the infall velocity based on the Zeldovich ap- 
proximation. The model was normalized to have 
the same infall velocity at the distance from the 
LG to the center of the Virgo cluster as found 
in the simulation. This model fares much bet- 
ter than the spherical Virgoccntric infall model of 
Huchra (1988): it makes good fit for all veloci- 
ties at all radii between the LG and the Virgo. 
Nevertheless, it failed because it required a mass 
between the LG and the Virgo that is about twice 
smaller than the corresponding mass in the sim- 
ulation (compared to the factor of three overpre- 
diction of mass by the Virgoccntric infall model) . 
Moreover, the velocity flow on the other side of the 
Virgo cluster at r > 15/i^^Mpc is not reproduced 
by this model even qualitatively. In fact, the veloc- 
ity flow around the Virgo cluster is clearly asym- 
metric and cannot be well described by a model 
with any kind of symmetry with respect to the 
Virgo cluster. This effort illustrates how difficult 
it is to construct a reasonable model for the veloc- 
ity fiow in the Local Supcrcluster. 

6. Discussion and Conclusions 

In the previous section we presented results of 
numerical simulations with initial conditions con- 
structed using large-scale constraints from the ob- 
served density field (deduced using the MARK III 
peculiar velocity survey) and a random realiza- 
tion of the density field for the popular ACDM 
model of structure formation. These initial con- 
ditions lead to the formation of large-scale (> 
5/i~"'^Mpc) structures that match the observed 
structures within ~ 100/i~^Mpc around the Local 
Group well, while the structures on smaller scales 
are not constrained and represent a random real- 
ization of the ACDM primordial power spectrum. 

We show that such constrained simulations re- 
produce the most prominent nearby large-scale 
structures, such as the Coma cluster, the Great 
Attractor, the Perseus-Pisces supercluster, the 
LSC, and the Local Void. The locations of the 
structures varies within the ensemble of the con- 
strained simulations (with different random real- 
izations of the ACDM spectrum) by a few mega- 
parsecs, but the overall qualitative large-scale 
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Table 1 
Peculiar velocity dispersions in the vicinity of the LG 



Distance from the LG 


K-adial 


V3D 


^halos 




conditions 


h-^Mpc 


km s^^ 


km s^^ 








5 


60 ±16 


129 


14 


Np> 


10 & Fcirc > 55 km s-i 


5 


58 ±16 


130 


13 


Np> 


15 


5 


57 ±16 


132 


12 


Np> 


20 & Fcirc > 60 km s-i 


5 


48 ±20 


136 


6 


Np> 


25 & Fcirc > 90 km s'^ 


6 


66 ±14 


139 


24 


N^> 


15 


6 


67 ±14 


141 


17 


N^> 


25 


8 


84 ±13 


165 


43 


Nj,> 


15 


8 


77 ±14 


181 


29 


Np> 


25 


10 


111 ±13 


222 


75 


Np> 


15 


10 


110 ±15 


222 


55 


Np> 


25 



morphology is reproduced robustly. The differ- 
ences in structure location from realization to 
realization imply that the constrained simula- 
tions should not be expected to reproduce the 
observed structures accurately in the highly non- 
linear regime. Instead, they should be considered 
as a tool for generating realizations of the primor- 
dial density field for specific purposes and compar- 
isons with observations in which it is crucial that 
the cosmic variance is greatly reduced. In other 
words, constrained simulations are well suited for 
studying the formation of particular objects, den- 
sity field configurations, and environments. The 
outputs of constrained simulations can be used 
for generation of mock galaxy catalogs and test- 
ing empirical models (e.g., the peculiar velocity 
field models). 

In this paper we used the constrained simu- 
lations (described in detail in § 3) to study the 
morphology of the density distribution and the 
peculiar velocity field in the ^ 30/i^^Mpc region 
around the Virgo cluster (which we call the LSC 
region in this study). In particular, we focused 
on the question of the relative "coldness" of the 
peculiar velocity field in the immediate vicinity of 
the LG. We find that the velocity field in the LSC 
region is quite complex. Matter and halos in un- 
derdense regions, such as the Local Void, undergo 
coherent expansion as predicted by analytic mod- 



els of void evolution (Hoffman & Shaham 1982; 
Bertschinger 1985). Near large-scale filaments and 
sheets the peculiar velocities are nearly perpen- 
dicular to the longest dimension of the closest 
structure, the pattern predicted by the Zel'dovich 
(1970) collapse solution. 

The main body of the LSC (the - 20ft."iMpc 
filament centered on the Virgo cluster) in the sim- 
ulation contains most of the mass in the LSC re- 
gion (« 7.5 X lO^^h'^Mpc), while the LG is lo- 
cated in an inconspicuous filament bordering the 
Local Void. The virial mass of the Virgo clus- 
ter is « 1 X lO^^h^^Mpc or less than 15% of the 
LSC mass and therefore the cluster does not dom- 
inate dynamics around the LSC. Correspondingly, 
a spherically symmetric Virgocentric infall model 
provides a very poor description of the peculiar 
velocity field which shows coherent infall in the 
overall direction of the extended mass concentra- 
tion of the LSC rather than towards the Virgo 
cluster. In particular, the model overestimates the 
mass within 10/i~^Mpc around the Virgo cluster 
in our simulation by a factor of three. The value 
of matter density Qq inferred using this model 
is thus grossly inaccurate. This inaccuracy was 
pointed out in previous tests of this model in the 
LG like environments in cosmological simulations 
(Cen 1994; Governato et al. 1997). 

The matter within and around the LG coun- 
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R(Mpc/h) 

Fig. 9. — Bottom panel: Peculiar radial velocities 
of DM particles (points) along the line-of-sight from 
the LG through the center of the Virgo cluster. All 
particles within the radius of 1.5/i^^Mpc of the line- 
of-sight are plotted. Solid circles with error bars show 
average and the rms deviation of the radial velocities 
as a function of radius; the rms radial velocity disper- 
sion is shown separately in the top panel. The Virgo 
Cluster is at the distance of 14/i~^Mpc from the LG. 
The particle velocities are computed in the restframe 
of the Virgo cluster. 



terpart in our simulation participates in the in- 
fall onto the LSC (roughly in the direction of the 
Virgo cluster) with a velocity of ^ 250 km s"-'^. 
Although the observational measurements of the 
LG infall velocity span a wide range of values, the 
average value, ^ 250±50kms^^ (see, e.g., Huchra 
1988, and references therein), is in good agree- 
ment with that found in our simulation. The local 
overdensity of dark matter around the simulated 
LG is w 5.5 within 1.5/i^^Mpc and is about zero 
within 5/i~^Mpc. This can be compared to the 
overdensity of galaxies of ^ 0.25 within 5/i~^Mpc 
in the IRAS survey (Schlegel et al. 1994). Given 
the small statistics and observational uncertainties 
this is a fair agreement: both the simulated and 
real LG appear to reside in only slightly overdense 
environments. 

Although the LG is located in a filament, this 
filament is of relatively small mass and does not 
perturb the large-scale infall flow signiflcantly. 



Therefore, the flow in the immediate vicinity 
of the LG is rather smooth. This explains the 
relative "coldness" of the local peculiar veloc- 
ity fleld. The radial velocity dispersion within 
w 7/i^^Mpc around the LG in the simulation is 
^ 40 — 70 km s^^, which is in good agreement with 
observations (e.g., Sandage 1986; Giraud 1986; 
Schlegel et al. 1994; Karachentsev & Makarov 
1996; Ekholm et al. 2001, and references therein). 
The velocity dispersion increases to « 110 km s~^ 
at larger distances due to larger peculiar velocities 
within and around the main body of the LSC. 

The low value of the velocity dispersion we find 
in our simulation is in contrast with considerably 
higher values found for the simulated LG like sys- 
tems in previous studies (Schlegel et al. 1994; Gov- 
ernato et al. 1997). We are not certain what 
caused the differences. We think that some of the 
differences are due to the different cosmological 
model. Very likely this is so in the case of the 
SCDM model as opposed to ACDM model stud- 
ied in this paper. In case of the open CDM model 
studied by Governato et al. (1997), which has the 
same density of matter as in our ACDM model, 
the differences may also be partially attributed to 
a more accurate representation of the large-scale 
environment of the LG in our simulation. In any 
case, it is comforting that the "cold" local veloc- 
ity field is reproduced within the current favored 
ACDM cosmological model and there is no need 
for exotic explanations (e.g., Chernin et al. 2001). 

To conclude, we think that numerical exper- 
iments conducted under controlled conditions of 
constrained realizations should help to shed light 
on some of the unanswered questions in the theory 
of structure and galaxy formation. Some of these 
questions arise from the unique data in the imme- 
diate neighborhood of the LG. The observational 
constraints on the primordial density field can be 
strengthened with the improving size and accuracy 
of the peculiar velocity surveys. Indeed, even in 
the immediate future one can conceivably perform 
more accurate constrained simulations by combin- 
ing the MARK III survey used in this study and 
the Surface Brightness Fluctuations (SBF; Tonry 
et al. 2001) peculiar velocity survey, which has 
more accurate velocities for nearby galaxies. 

The question of the coldness of the local pe- 
culiar velocity field addressed in this paper is 
but a single example of the problems that can 
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be addressed by constrained simulations. Indeed, 
in a related paper we extend the present work 
by using gas-dynamics constrained simulations to 
study properties of the intergalactic medium in 
the LSC region and formation of the Virgo cluster 
(Kravtsov et al. 2001). Wc also plan to use higher 
mass and force resolution constrained simulations 
to study the distribution of dwarf galaxies in the 
LSC and nearby voids, formation of the LG and 
other nearby groups and clusters, etc. Overall, we 
find that the constrained simulations provide an 
optimal tool for studying a wide range of dynami- 
cal problems concerning our "local" neighborhood 
within the large-scale structure of the Universe. 
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